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ABSTRACT 


STICAP (Stiff Circuit Analysis Program) is a FORTRAN IV, Version 
2.3, computer program written for the CDC-6400-6600 computer series 
and SCOPE 3.0 operating svstem. It provides the circuit analyst a 
tool for automatically computing the transient responses and fre- 
quency responses of large linear time invariant networks, both stiff 
and non-stiff. The circuit description and user's program input 
language is engineer-oriented, making simple the task of using the 
program. 

Three volumes of documentation are available for the STICAP 
program; a theory manual, a user's manual, and a system's programmers 
manual. Volume I describes the engineering theories underlying 
STICAP and gives further references to the literature* Volume II, 
the user's manual r explains user interaction with the program and 
gives results of typical circuit design applications. Volume III 
depicts the program structure from a system's programmers viewpoint 
and contains flow charts and other software documentation. 
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CLOSER I 


GENERAL PRCGRAII DESCRIPTION 


1 . 0 INTRODUCTION 

The program STICAP - ( Stiff Circuit Analysis Program - was 
developed by personnel of the School of Engineering, Old Dominion 
University, Norfolk, Virginia, 1970-1971, under contract NAS1-9434-25. 
This program package represents the merging into one diversified com- 
puter aided network design program tne capabilities of the existing 
programs CORNAP, 1 for linear circuit analysis; Gear’s ALGORITKK 407 - 

O 

DIFSUB,' 1, for numerical integration of stiff ordinary differential 
equations? and a somewhat specialised matrix solution technique for 
obtaining time domain circuit response. 

The composite program thus consists of three separate component 
programs, or modes of operation, each with some advantages over the 
others, in different circumstances. The COBITAP mode consists of the 
circuit analysis programs and capabilities of the original program 
CORNAP. In the Gear and matrix modes the circuit translation routines 
of the program CORNAP are employed to obtain the state variable dif- 
ferential equations of the circuit, but different techniques for 
solving these equations are used. The functions and limitations of 
the various inodes are described in the sequel. 

The program STICAP is written m the FORTRAN IV, version 2.3 
language. It is machine compatible with the CDC 6400-6600 computer 
series and runs under the SCOPE 3.0 operating system. It is segmented 

^Developed by Dr, Christopher Pottle, Cornell University, Ithaca, N.Y. 
^Developed by Dr. C. V 1 . Gear, University of Illinois, Urbana, Illinois 
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xn overlays of 70K or less using the SCOPE OVERLAY capabxlity. All 
I/O xs accomplished using standard I/O files . The I/C files are 
equivalenced so that File' 5 is used for input and File 6 for output. 
Mo other files are used by this program. 

1.1 PROGRAM FUNCTIONS AMD CAPACITY 

This program has the capability of obtaining, at tne option of 
the user, certain combinations of the following quantities; state 
variable equations, transfer functions, frequency and time responses, 
of an n-port linear active time invariant network. 

The starting point for the program’s analysis is a user oriented 
circuit description stated m rerms of circuit branch elements and 
circuit nodes. The largest network configuration of these elements 
accepted by the program may be determined as follows- Let E be the 
number of energy storage elements ; I the number of inputs ; IT the num- 
ber of outputs; R tne number of resistors; and C the number of con- 
trolled sources present m the network. The maximum number of ele- 
ments of each type are governed by the constraints 

E + I < 30, 

E + TT £ 30 , 

R + C < 28. 

In the Gear and Matrix modes the additional constraints 

1 < 1° , M < 10 

are imposed. 

1.2 PROGRAM SELECT OPTIONS 

The user may select one of tne following mutually exclusive 
modes of operation; the CO REAP mode, the GEAR mode, or the matrix 



mode. The first of these, the CORHAP mode, embodies the network 
analysis- capabilities of the original nrograin OORHAP , with choice 
suppression of certain print features, such as the printing of state 
equations and transfer functions. This program mode translates a 
circuit description from user language in terms of circuit nodes and 
branch elements to a mathematical description m terms of the state 
variable differential equations and algebraic state-output equations 
of the circuit. The option is provided for subsequent calculation 
of transfer functions, zeroes of transmission, and frequency or time 
response of the circuit. 

The Gear mode may be used to perform time response calculations 
only. Here the circuit equations are generated by the CORNAP sub- 
routines, and either stiffly stable implicit linear multistep methods 
or the non-stiff Adam’s integration techniques may be selected for 
numerical integration of the state equations. In this mode a maximum 
of ten independent sources nay be simultaneously used to drive the 
network, and up to ten simultaneous outpurs may be requested. The 
full power of tne F0PTRA2-I language may be employed to describe the 
mathematical equations governing the behavior of the independent 
sources , or the user may wish to write his own program for input of 
sampled data. 

Finally, the matrix mode may also be used for performing time 
response calculations, employing the circuit equations generated by 
CORNAP subroutines . The solution of these equations is obtained by 
means of a matrix technique which avoids a numerical integration. 

The technique is computationally rapid, but it is applicable only m 
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the case of linear time invariant systems whose eigenvalues are not 
closely- grouped, and ”hich are forced by sinusoidal, cosinusoidal, 
or step function inputs. Only a limited number of such inputs are 
allowed. The circuit may ae driven by a maximum of ten simultaneous 
independent sources , and a maximum of ten outputs may be requested. 

1.2 WETUORK ACCEPTIBILITY 

This program package will perform the complete analysis of any 
lumped, linear time invariant network, whether stiff or non~stiff. 

The elements making up the network nav be of the following types: 

1) ordinary two- terminal passive circuit elements - resistance, 
inductance, and capacitance; 

2) mutual inductance and capacitance; 

3) the four two-terminal controlled sources (voltage/current 
controlled voltage/current sources) . 

Two port active and nonreciprocal elements such as negative impedance 
converters, ideal transformers and gyrators can be made up of rile one 
port elements described above. Inputs are defined by attaching in- 
dependent voltage and current sources to the network. Unity coupled 
trails formers (or even n-port inductors with a semidefimte inductance 
matrix) can be handled by the procedure, as can all resistive networks 

1,4 SOURCE DERIVATIVES 

Source derivatives will occur m the state variable equations 
describing a circuit whenever a voltage source is connected m a loop 
containing only capacitors and other voltage sources, or a current 
source is connected m a cut set containing only inductors and other 



current sources. Although the state variable equations describing 
passive circuits nay contain only the first derivative of circuit 
inputs, the equations describing active circuits iray contain any 
number of input derivatives. An active circuit whose state variable 
equation contains a second derivative is shos^n m Figure 1. The 
state variable and output equations for tins circuit are 


dv c3 

“dt" 


“•5V c3 + 


_ cl 2 V, 

5 c 

dfc^ 


v 0 ~ v c3 


A. 

dt^ 


I The CORHAP mode will nor compute jnne domain responses for 
circuits with xnpur derivatives in the state output equations . 

The GFA.R and IL'TPIR . odes ello^ only a first derivative of tue 
circuit inputs m the state variable equations as long as no input 
derivatives occur m the output equations, in this case the state 
variable equations and the output eouations vrould be of the form 


X * AX + Bu + B x u 
Y = CZ + Du 

A second choice of state variables 

a = X + B^u 

would then be made transforming the original state variable and out- 


put equation into 




Figure 1 
L 



Figure 3 
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where no input derivatives occur. For example, consider the circuit 
m Figure 2. The state variable eauations for this circuit are 



TThether an input derivative occurs in the output equations depends 
upon the choice of the outputs. If v c3 were chosen as the output, 
the output equation would be 


V 


c2 


= EO-1] 


fx L l 

y 


+ v. 


For the GEAR and the J5A.TRIX mode, these equations would then be trans 
font>ed into the ecuivalent form 


V 
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If, however, i c were chosen as the circuit outout, the output equa- 


tion would be 

^■cl = ([~7f'°) 


+ 

L v cij 


_ dV 3 

4 i 

dt 


and no solution would be attempted, bv any of the three solution modes 


Care should be taken when describing the inputs to a circuit 
with input derivatives. Due to the fact that inductor currents and 



8 


capacitor voltages cannot change instantaneously , inputs to these 
types of cxrcuxts must not be allowed to have ;jump discontinuities. 
Jioreover, the xnxtxal conditions of the circuit nust be chosen such 
that the cxrcuxt at time t = 0 + obeys "xrchhoff * s voltage and current 
laws. For example, xf the cxrcuxt xn Figure 2 were to be driven with 
the voltage source 

V x = cos (t) 

the initial capacxtor voltages must be cnosen such that 

Vcl(O) + V c2 {0) = V x (0) « 1. | 

In that the initial conditions of a circuit are chosen to be sero 
for impulse and step calculations, neither the GEAR nor the IIATFIX 
mode will calculate the step or the impulse response of a circuit 
whose solution equations contain input derivatives. 

The source derivatives can be eliminated from the solution equa- 
tions by including the internal resistance of the voltage sources 
and the parallel conductances of the current sources. Thus, if the 
internal resistance Rg of the voltage source in Figure 2 were in- 
cluded m the circuit description as sho^m m Figure 3 , the state 
variable and output eauations would become 




CHAPTER II 


CIRCUIT DESCRIPTION MD HODS SELECTION 


2.0 OVERVIEW OF CARD INPUT DECK SETUP 

In this chapter the manner m which the user describes his 
circuit to the program STICAP X7ill be discussed. As a means of 
introduction an illustration of the overall deck setup of the cards 
which must be prepared is indicated by Figure 1. 




END CARD 


Data Cards for Sam- 
pled Inputs; CORNAP 
mode only 


Control Cards for 
the node Selected 


Hode 

Select 

Card 


r~ 


Scaling 

Card 


Outputs 

Description 
Cards Group 


Elements Descriptior 
Card Group 



Figure 1 - Overview of Deck Setup 
The contents of the title card are printed out verbatim at the 
head of each section of the output, and serve as a means of identi- 
fication to the user. As such, the user has complete freedom m 
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specifying tne contents of Jus curd. Tne elements description card 
group is usel co specify the interconnections neeween nodes and 
branches of tne circuit. L'ne oxitpucs card group is used co specify 
the circuit currents anc. volteigcs cesirad as outputs w tne user. 

Tne scaling card specifies the n aimer m wmch tno circuit is to be 
jcalec for computational purposes. The mode selecc card determines 
v.’nicn of tne Gear, CORNAP or Matrix modes is to be selected for 
analysis. A different group of control cards is needed, depending 
upon vmch mode is selected. These cards specify tne output options 
selected and contain information needed ny the numerical integration 
and analysis routines. In che COKHAP irode only, a series of card 
groups consisting of a control card followed by sampled input data 
cards may sometimes appear. An end card is always present, regard" 
less of mode. 

2.1 ELEMENTS DESCRIPTION CARD GROUP 

The first card m tne group is a header card 
* SLEuEMTS 

containing an asterish m column one. Tne word ‘ELSiiJMTS" may 
appear anywhere on the card (data fields whose lengtn and starting 
position on tne card may be arbitrarily selected by tne user are 
said to be free form oata) . 

All elements description cares are free form. One card must 
ns prepared for each circuit element. Card format for passive 
elements or independent sources is 

NAME HI J2 VALUE, 

ORIGINAL PAGE IS 
OP POOR QUALITY? 
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ana for independent sources xs 

UAME Hi M2 V/01EJHTAME2. 

HAMS, Ml, IJ2, and VALUE must be separated oy at least one blank. 
The asterisk between VALUE and HAMD2 need noc be present. 

ILuIE is the user ! s name for the circuit element. T_ns name 
nay be a maximum of four characters long. The first letter of tne 
element name is used to specify the element type as indicated s 
first Letter Element Type 

V Voltage source 

I Current source 

R Resistor 


a 


Inductor 


C Capacitor 

K Coefficient of Coupling 

M Mutual inductance or Capacitance 


For non-mutual elements Hi and W2 are the node numbers (two 
digit integers) of the nodes between wmch the element is connected. 
The circuit r.oaes snould be ordered compactly from zero (00) , altnough 
the failure to do so is non- fatal. Hode 00 should be the reference 
node, or ground. The maximum nu’iher of nodes is 54. HI is given 
a positive reference with respect to H2“ 



For mutual elements HI anu IJ2 are the element names of the two element 
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mvolved. 

For dependent sources NAME 2 is tlxe name of the controllxng ele- 
ment, i7i tn a V or an I prefixed as the firsc letter of the name, 
to indicate respectively tne occurrences of a voltage controlled 
or current controlled element, 

A feature not normally employed by the average user is the 
following- By placing tne word TRUE 1 after the description of the 
element, as m 

NAME N1 N2 VALUE TREE, 

a capacitor may be forced into the proper tree. An inductor may 
similarly oe forced into the cotree oy placing the word " COTREE " 
there, i.e., 

ITAUE 151 LT2 VALUE COTREE. 

VALUE is tne value of tne element or strength of the source. 
Negative values are permissible. Zero values are also allowed for 
inductors and capacitors and may be used to define fictitious branenes 
for oueput or control purposes. No entry need be given for indepen- 
dent sources. VALUE may be any mceger, decimal, or floating point 
number with a maximum of 15 digits m the mantissa, and where mag- 
nitude is m the range 1Q“ 230 to 10 +290 . 

2.2 OUTPUTS DESCRIPTION CARDS CROUP 

Tne element currents and voltages selected by the user as circuit 
outputs are indicated oy tnxs card group. Tne neacer card, free 
form excepc for the asterisk m column one, has the format; 




* OUTPUTS 
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The output cards, ordered m anv fashion, are of oae of the forms 
VNAME 

or INANE , 

depending upon whether the voltage across or current through the 
element with this element name is desired. Here NAME is the name 
of a circuit element given in the elements description card group. 

One card must be prepared for each desired output. 

2.3 SCALING CARDS GROUP 

The scaling cards, in free format, are preceded by a header 
card with an asterisk m column one 
* SCALING 

followed by one or both of the cards 
FREQUENCY = VALUE! 
and IMPEDANCE = VALUE 2. 

VALUE 1 and VALUE 2 represent respectively the frequency (rad/sec) 
and impedance level (ohms) about which the network is designed to 
operate. These numbers, ideally, are used as scale factors to scale 
the network to operate around 1 rad/sec and a 1 ohn impedance level. 

A. scaling factor several magnitudes away from its proper value can be 
computationally critical? hence some attention should be paid to 
determining scale factors which cause scaled element values to brac- 
ket the value umtyt The scaling values nay be integer, decimal, or 
floating point, as specified in section (2.1) , 

The scaling cards group ray be omitted, but if so the user 
should be careful to use a consistent set of units for the circuit. 
Most standard texts on network theory include a section on circuit 


*See Appendix I for further discussion of scaling. 



scaling, and the concept should he understood by the user before 
attempting to design with STICAP. 

If the scaling option is chosen, the Gear and >!atnx mode output 
remains scaled by the same factors. In the CORPAP mode the choice 
of unsealing the output is present. 

2.4 HOPE SELECT CARD 

Regardless of which mode is to be selected, the user's deck set- 
up is the same as is indicated by Figure 1, up to and including the 
scaling cards. Thereafter it changes with the individual mode. In 
this section the mode select cards and the control cards for each 
mode will be described. 

The mode select cards are of the general form 

*flODE NAI5E, OPTION, OPTION, OPTION. 

The asterisk necessarily appears m column 1, followed by the mode 
name of the mode selected, and the options chosen. All data fields 
other than the asterisk may be free form, but must: be delimited by 
commas. The omission of any option indicator results in the omissio] 
of that particular option. The options need not occur m any speci- 
fxc order on the mode select card. 

The mode select cards for each mode, with all possible options 
present, are given below. 


IIODE HAtZ 

OPTION INDICATORS 

CORNAP MODE 

STATE EQUATIONS, TRANSFER FUNCTIONS 


STATE EQUATIONS, TRANSFER FUNCTIONS 

IIATRIX IIODE 

STATE EQUATIONS, TRANSFER FUNCTIONS, 
SOLUTION EQUATIONS 
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In all three nodes the matrices of the state-input and state output 
equations 

X* =s 'Ax + Bu + ciu* + B 2 U 1 + . . . 
y = Cx + Du + E^u* + D 2 u ,f + . . . 

nay be printed, using the option indicator "STATE EQUATIONS." How- 
ever, if a higher order derivative than the first occurs m the 
state-input equations, or if any derivative at all occurs m the 
outout ecruations, numerical integration of these equations cannot 
be performed, since no allowance is made for user input of source 
derivatives „ 

If the transfer functions indicator appears on the node select 
card, the poles, zeroes, and gam constants of the transfer functions 
of each input-output pair are printed. In the natnx mode the net- 
work ecruations are presumedly solvable m closed form, and the exact 
solution equations may be printed by means of the third option indi- 
cator. This feature prints the solution equations governing the time 
domain behavior of the state variables . The combined output of these 
ecruations together with the matrices of the state output equations 
allows the user a complete closed form description off the solution. 

Example 1 

The following cards effect the same result, the choice of the 
Gear node with its full options : 

* GEAR HODS, TRAITSF3R FUITCTIOHS , STATE EQUATIONS 
*GEAR, STATE, TRAP'S FER 
*G, STA, TRA H 
*G , TRA , STA 
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The program reads only the first letter of the mode name and the 
first 3 letters of each option indicator. Errors in the mode name 
are fatal? those m options selection are merely indicated by error 
messages, with the indeterminate option omitted. 

Example 2 - Fatal error m mode name. 

*HEAR 

Example 3 - Error diagnostics; liATRIF option chosen? contrary to 
desires of the user; print request ignored. 

^IlEAR TRANSFER (Comma omitted) 

*HORNAP ,T RANSFDR (Incorrect option indicator) 
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chap ' ms r hi 


CONTROL CLODS GEAR hODE 


3.0 GENERAL NODE DESCRIPTION 

The prmary purpose of the ^eer program mode is the obtaining 
of tire domain circuit response for stiff circuits , using stiffly 
stable implicit linear nultistep methods for perforiring the nur>crxcal 
integration. The response of non-stiff circuits may be obtained 
eciually as ^ell, since the option of choosxnc the Adair's linear mul 
tistep methods, suitable for non-scif f integration, is available. 
Since little extra effort is required to do so, the op-cxonal capabil- 
ity of obtaining the poles, zeroes, and gam constants for tne trans- 
fer functions of each input output oair m the network is included. 
Tne tire domain response of a circuit with up to ten simultaneous 
inputs or outputs nav be obtained. The inputs nay be described using 
the full power of the FORTRAN language, or the user may wish to write 
his o'm routine for inputting sampled data, as indicated in the 
description of the user routine USEFCM.* Thus the Gear mode is the 
most powerful general purpose analysis component. Its one possible 
drawback, as opposed to bhe COKNAP mode, is that the input waveforms 
cannot be changed without reprocessing the entire circuit, anc it 
cannot be used, for frequency response calculations. 

The control cards for this node T ' 7 ill now be described {see 

Figure II) . 


* Sub routine USEFCN is a user supplied routine which specifies the 
independent sources of the circuit, discussed at a later point. 
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End Card 


1 Hun Controls 
Cards Group 


Initial Conditions] 
Cards Groum 


Source Order 
Cards Group 


Figure XI - Control Cards ? Gear I'ode 
These control cards may be classed m the categories s source order 
cards, initial conditions cards, and run control cards. It is not 
essential that the ordering of the three card groups within the deck 
setup be as indicated by Figure II j however, this ordering should 
process fastest. Cards comnosmg each individual group will now be 
described. 

3.1 SOURCE ORDER CARDS GROUP 

(a) Header card - The first card m this group contains an 
asterisk m column 1, followed by the words ’SOURCE ORDER- 
IMG,” m free form. 

(b) Source name cards - These cards must be ordered ±n the 

seguence that the values of the independent source functions 
U(l), U(2) .. r U(3f) are commuted by the user-supplied 

routine USEFCH. The name of the corresponding independent 
source aopears any /here on the card. This name must be 
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the same as the naire given this source on the circuit 
description cards. 

Ex 



The number of source names must be the same as the number 
of independent sources defined when describing the circuit, 
else the program is terminated. If the number of names 
does not agree with the number of functions defined m 
USEFCN as sources , the program is also term mated. If an 
impulse response is to be computed, only one independent 
source is involved; hence, only one source should be speci 
fied when describing the circuit. 

The source names and the order these source values 
are computed is necessary m order to set up the proper 
correspondence between components of the source vector as 
computed bv the user, and the order in which these sources 
appear ordered by the circuit translation program. 

If only one source is defined in the circuit descrip- 
tion, this card group need not be present. 

3.2 INITIAL CONDITIONS CARD GROUP 

(a) Header card - The first card m this group contains an 
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asterisk in column 1 , followed bv the sort's Initial Con- 
ditions, 1 m free form; i„e. no specific starting column 
for either word. 

(b) Value cards - These cards may be ordered m any fashion and 
are also free form, Hach card contains the name of a cor 1 
ponent of t^e state vector of the circuit, as well as its 
value, m the form 

VATU = Value. 

Here I7AHH is the same as the name specified on the cir 
cuit description cards, and ’Value has one of the formats 

i 

a) integer 

b) decimal 

c) exponential - exponent field of length 3. 

The components of the state vector are the capacitor vol 
tages and inductor currents of the circuit. 

Examples - 

Cl = .003 
1.2 = 2.S + 003 
C2 ~ 3E - 0 03 
C5 = 6 

If this card group is not present, it is ass^^med that the 
initial state is sero. Any state variable not assigned an 
initial value is assumed to have an initial value of sero. 

3.3 HUH CONTROLS CARD 

This card group supplies data needed by Gear's program, and 
allots selection of certain options. 


ORlGIHMiPAfElS 

OF POOB 


Individual cards m the 
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group are as follows; 

(a} Header card - The first card xn the group, whxch contains 
an asterisk xn column 1, follwed by the vzords '’RUN CON- 
TROLS," in free fom. 

(b) Integration control cards - These cards my be ordered in. 
any fashion and are m free form. One card is required for 
each of the following items , if the item is selected as an 
option by the user; 

1. Initial time - ^his is the lower limit of integracion 
and the time at which the initial conditions are 
measured. This card need not be given if the initial 
time is zero. It will be of the form; 

INITIAL TIME = A 

where A is a floating point number, 

2. Impulse or step response - The impulse or step response 
of the circuit is found if one of the following two 
cards is present; 

STEP RESPONSE 
IMPULSE RESPONSE 

Only one of these cards nay be given. Note that 
USEFCil is not to be supplied by the user if the im- 
pulse or step response is to be calculated, 

3 „ Print starting time - This card gives the value of time 
at which output printing is to begin. The card will 
be of the form; 


PRINT START « A 



SSKfiSS 
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where A is a floating point number. If this card is 
not given, printing s will begin at the initial tjne, 

4. Upper integration limit ~ The upper limit for the 
integration nay be given by specifying the stop time 
or bv specifying the number of time points printed* 

This card will be one of the t"’0 forms z 

STOP TIMS = A 
POINTS PRINTED - 17 

where A is a floating point number and n is an integer 
number. If a stop time is given , interpolation is 
used to determine the values of the outputs an the 
stop time. If neither card is given, 100 time points 
will be printed. 

5. Error controls - The value of EPS controls the accuracy 
of Gear's integration routines. The Euclidean norm of 
a vector whose Ith component is the single step error 
of the Ith state variable divided by the maximum value 
of the Ith state variable, must be less than tins 
value. Tins card trill be of the form 

EPS = A 

where A is a floating point number. If this card is 

- 4 

not given, EPS is assumed to be 10 

6. Integration method - An integration method suitable 
for stiff systems is normally used. A predictor-cor- 
rector Adam's integration method, however, may be 
chosen by including the following cards 

ADA.4S INTEGRATION 
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i 

7. Print density - The number of integration steps between 
printings is given by this card and is of the form 

OUTPUT DENSITY = N 

where IT is an integer. If not given, printings will 
take place after every third integration srep. 

8. Step size controls - The size of tne integration step 
is controlled bv the integration routine," however, the 
maximum step size, minimum step size, and initial step 
size may be specified by the user. These cards will 
be as follows * 


HKAX = A 
Hi 1IU ** A 
HiniT = A 

where A is a floating point number. If HI1AX is not 
given, it is set to one tenth of the stop time minus 
the initial tine, if the stop tine was specifsed, and 
is unbounded otherwise. If HI HR is not given, in is 
ser to 2 ero„ If HIHIT is nor given, it is set to 10 


3 . 4 END CARD 

The last card m the Gear input is an END card. The format is: 

* END 

with the asterisk in column 1. 

3.5 USER SUPPLIED INPUT ROUTINE 

The routine USEFCN is a subroutine which must be supplied by the 
user. The function of this routine is the following s Given a 
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specific time at which their values are needed by the integration 
routine, compute values U{1), I = 1,2 , . „ .N,II<10 of the various m~ 
dependent sources and store them in an output source vector U„ 

Tms feature of STICAP supplies to the user two alternative 
capabilities ; First, the independent sources may be defined m 
equation form, using the full power of the FORTRAN language. Second, 
values of the independent sources may be obtained by interpolation 
of data samples. This data could be supplied in tabular form, or be 
read m block by block at execution time. 

Tne program STICAP has incorporated within it a skeleton USEFCM. 
This routine may be completed by addition of user supplied FORTRAN 
equation statements defining the independent source equations ; or it 
may be replaced by a user supplied routine which at specified times 
computes by some means, such as indicated, the vector tJ of source 
values. Assuming STICAP resides on a data cell, this routine could 
be caused to replace the resident USEFCN routine prior to execution, 
by means of a CUTOUT card, or whatever UPDATE facilities are available. 

It is essential that the source values m the USEFCII output 
vector be ordered xn the same sequence m which the names of the 
independent sources occur on the user's SOURCE ORDERING cards. "Te 
emphasize that errors m this ordering cannot be detected by STICAP, 
and result in a grossly misleading circuit analysis. 

Sources defined by FORTRAN equations 

If there are N independent sources , the user supplies N state- 
ments of the form 

U(l) « F1(T) 

U (2) = F2(T) 
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Here T is a dummy argument whose value , when sunplied by STICAP , 
specifies the time at which a vector of source values is required. 

The FI (T) , I - 1,2,..., ft, are FORTRAN equation statements defining 
the sources as a function of the dummy time variable T. The ordering 
of the values of U(I) is identical to the order m which the names 
of the sources occur on the source order cards. 

The skeleton USEFCM contains the soatemen-cs below; 

SUBROUTINE USEFCft (T,U) 

DIMENSION U(10) 

COI S ION /RUN/SKI 13 (13) , INPUT 
IF (IftPUT-EO-1) GO TO 10 
U(l) = 1.0 
GO TO 20 

10 U(l) = 0.0 

20 CONTINUE 
RETURN 
END 

The FORTRAN equation defining statements produced by the user 
are to replace all statements starting with the COITION statement and 
ending with statement number 20. These statements may be replacer3 
using whatever UPDATE facilities are convenient. 

Source Values Obtained by Interpolation 

The circuit designer may wish to supply an interpolation routine 
which obtains the independent source values from tabular data de~ 
fined within itself, or from data samoles read in at execution tine. 
If so, the onus is on the user to determine the necessary data den- 
sity and precision of the interpolation scheme recmired in order to 
provide accurate intermediate values. Here it should be kept m 
mind that data values needed by Gear's integration scheme are not 
uniformly spaced, nor do they necessarily occur m a time sequence 
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which is ordered m terms of increasing time. Hence data points as 
much as- one plot point (print increment) behind the current time 
point may next be recmired. 

If this USEFCN option is chosen, the user*s program size must 
be compatible with overlay sizes of the mam program. The user 
written routine may be inserted m the program replacing the above 
listed USEFCN using whatever UPDATE facilities are convenient to the 
user. Any data cards read by this routine should follov/ the END 

card in the deck setup indicated by Figure I, Chapter II. If the 

l 

data is read from tape, the tape cannot be defined as logical unit 
"5 or #6. The program card m the mam overlay should be modified 
to include this tape number. 
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CONTROL CARDS, MATRIX MODE 
4.0 GENERAL MODE DESCRIPTION 

The Matrix program made provides a rapid means for obtaining tine 
domain circuit response, for circuits with inputs restricted to the 
class of linear combinations of smusoidal, cosinusoidal, step and 
impulse functions. Stiffness of tae circuit does not affect the 
analysis, but rainier the grouping of eigenvalues of the system matrix. 
Such eigenvalues should not be too closely grouped together, or 
computational error may become a problem. The optional capability 
of obtaining tne poles, zeroes and gam constants for the transfer 
functions of each input-output pair m tne network is included. 

Input waveforms may not be changed without reprocessing the entire 
circuit, and no frequency response calculations may be performed. 

The control card deck setup for this mode is similar to that 
of the Gear mode (see Figure III) . It is not necessary that the 
ordering of the card groups conform to tne one given, however, this 
ordering should process most rapidly. 


END 

CARD 




Source 
Definitions 
Card Group 


'Run Controls 
Card Group 


f Initial 

Conditions 
Card Group 


Figure III - Control Cards; Matrix Mode 
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Tne composition of each card group will now be given . 

4.1 INITIAL CONDITIONS CARD GROUP 

The cards needed m this group are the same as those of the 
Gear mode initial condition's description (see section III. 3. 2). 

This card group is not required if all initial conditions are zero, 
as m impulse response calculations. 

4.2 RUN CONTROLS CARD GROUP 

Individual cards in this group are as follows: 

(a) Header Card - The first card iiji the group, which contains 
an asterisk in column one, followed by the words "RUN 
CONTROLS", in free form. 

(b) Run Controls Cards - These cards may be ordered in any 
fashion, and may be punched m free form. One card must be 
present for each of the following options desired by the 
user: 

1. Initial time - the lower limit of integration, and time 
at which initial conditions must be measured. Card 
format is 

INITIAL TIIIE » A 

where A is a floating point number. The initial time is 
assumed zero if this card is omitted. 

2. Response type - the impulse or step response (with a 
step of amplitude 1) is calculated if one of the fol- 
lowing format statements is present: 

STEP RESPONSE 


IMPULSE RESPONSE 
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If neither card xs present, a source definition card group 
- specifying the network inputs is mandatory. 

3. Print starting time ~ Contains the time value at which out- 
put printing is to begin. Card format is 

PRINT START = A 

where A is a floating point number, if this card is omitted, 
printing begins at the initial time. 

4. Upper limit of integration - Hay be specified by a stop 
time or number of points to be printed. Card format is one 
of the following: 

STOP TIME - A 

POINTS PRINTED = N 

1 

where A is floating point format and N is integer. If 
neither card is present, 100 time values of the outputs will 
be printed. 

5. Print density - The plot increment between print points is 
specified in the format 

PLOT INCREMENT = A 
where A is floating point format. 

4.3 SOURCE DEFINITIONS CARD GROUP 

This, card group consists of a header card followed by a source 
definition card for each driving source. The header card contains an 
asterisk in column one, followed by the words "SOURCE DEFINITIONS CARD 
GROUP", in free form. A maximum of ten independent driving sources is 
permitted. Each source must be specified by a card (or cards) punched 
free form m the following format 

NAME * F1(T) + F2(T) + .... + FI>1(T) . 
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Here NAME is the name of an independent source specified when des- 
cribing* the circuit. The right member of the equation statement is 
a sura of functions FX(T), I <20, where each FI(T) may have one of the 
forms 

A 

A* IMP 
A*SXN(B*T) 

A*SXN(B*T+C) 

A*COS (B*T) 

A*COS (B*T+C) 

where A, B, C are free format floating point numbers. 

The first form indicates a step function of amplitude specified 
by A, with jump at the initial time t c ? the second on impulse func- 
tion A6 (t-t c ) , the delta function of amplitude A. All other forms 
indicate sines and cosines of amplitude A, frequency f , 2‘tif : =B, and 
phase C. The asterisks indicating multiplication need not appear. 

A source description may be continued on additional cards by placing 
a dollar sign $ in column 1 of each continuation card, as long as 
each FX(T) is completely described on one card. No more than 20 FI(T) 
may compose one source; a sinusoid or cosinusoid with non-zero phase 
is considered as two functions, for purposes of counting. 

4 . 4 END CARD 

The end of the source definitions card group, and matrix input, 
is signalled by a card with an asterisk m column 1, of the form 

* END. 
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CONTROL CARDS; CORNAP MODE 

5.0 GENERAL MODE DESCRIPTION 

The CORNAP mode embodies the capabilities of the original pro- 

I 

gram CORNAP. This program mode may be used to obtain transfer func- 
tions, zeroes of transmission, and frequency or time response of 
the network. Assuming a network characterized by multiple indepen- 
dent driving sources and/or multiple output ports, the user may obtain, 
for each input-output pair, the preceding quantities. Furthermore, 
it is not necessary to reprocess the circuit description in order 
to obtain the outputs at a fixed port which occur when a different 
input port is used to drive the circuit. The same is true if it 
is desired to alter the wave form at a fixed port. Unfortunately , 
no capability is provided for computing a superposition of the out- 
puts at a single port, assuming it were desired to drive a network 
wx th several simultaneous inputs. However, both the Gear and Matrix 
inodes have this capability. For time domain analysis, step response, 
impulse response, or transient response with a sampled data driving 
input may be obtained. In all cases it is assumed that the initial 
state is zero. However, m both of the other modes, the initial 
state vector may be selected arbitrarily, m these instances for which 
such choice might be desirable. 

5.1 CONTROL CARDS; DATA CARDS 

The cards following the mode select card (see Fig. 1; Chapter 
II) needed to complete the user’s deck setup for the CORNAP mode will 
be described m this section. These cards are not free format; infor- 
mation appears in a specific data field in a specific form. These 
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cards control the calculation and printing of time and frequency 
response data. In reading the card format description given below, 
it would be advantageous to keep m mind the following items: in 

the circuit description for the CORN&P mode, the network may be 
described as having multiple input ports and multiple output ports; 
the circuit translation routines then yield the proper set of state 
and state output equations for such a network. However , in the actual 
computation of circuit response, the responses are available m terms 
of input-output pairs; i.e, , the response to a single input- single 
output linear time invariant system is computed. The C0RNA1? mode 
does not compute a summed response at a single output port for a cir- 
cuit excited by several simultaneous inputs. Responses may be 
obtained for all possible input-output pairs/ without reprocessing 
the circuit. 

The control cards are prepared in the format below. However, 
if only state equations and transfer functions are desired, the cards 
now to be described may be omitted. 

Any number of these cards may be present, in any order. If 
sampled data input is desired, one of these cards must be present 
for each individual desired input, with the data samples immediately 
following, as described below. 

Col. 1 

This column contains a character which defines the type of 

response desired: 

F - frequency response 


T - time response. 
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These columns contain eit-ier the namli of an xnout-defining 
independent source or are bianjc. Responses with this element 


as input are calculated if a name appears here? responses x^ith 
eacli and every independent source as input are generated if 
the field is lefu blank. The field must not be left blank. 


for sampled data time response. Similar columns for output are 
cols. 33-36. 


Cols. S-9 


Frequency Response 

These columns contain a tuo-digifc integer giving the number of 
decades of the frequency variable to be covered by a frequency 
response calculation if a logarithmic scale for frequency is 
chosen. Blanks m these columns indicate a linear scale is 
desired. 


Cols. 10-12 


A three-digit integer m these columns give the number of fre- 
quency or time points to be printed. For a logarithmic fre- 
quency scale, this number gives tne number of points per decade . 
Cols. 14-31 

The increment m frequency or time netv/een printed responses 
appears in these columns as a string of digits containing a 
decimal point, Ssrpcnentxal noration (I.5E -04 = 0.00015) may 
be used provided the exponent part is nght- 3 ustified m the 
field (cols. 2S-31) . 

Cols. 33-36 

Tnese columns contain either tne name of an element defined in 


*left justified 
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the outputs description (without the appended V or I) as deter- 
mining an output port or are blank. Responses with this element 
as output are calculated if a name appears here? responses at 
each and every defined output are generated if the field is 
left blank. 

Col. 30 

Any nonblank character m this column will cause response calcu- 
lations to remain scaled by the factors given previously. These 
scale factors will be used to "unnormalize" the calculations if 
this column is left blank. 

Col. 41 

Frequency response 

A nonblank character in this column will cause the frequency 
scale to be m radians/sec. A blank m this column indicates 
the frequency scale is to be m Hz . 

Tima response 

A blank xn this column indicates the impulse and step response 
of the network are to be calculated? a nonblank character causes 
an external sampled input signal to be used as input. At present 
the samples of these signals are entered 6 to a card, each 
occupying 12 columns (cols. 1-12, 13-24 ,61-72) , m the same 
way as the time increment between printouts was entered. As many 
input signals must be entered as inputs were defined in cols. 3-6 
The first sample of an input signal must begin a net*? card. These 
cards follow immediately the time response control card now tinder 
discussion. 

Cols. 44-61 
Frequency response 
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The first (lowest) frequency at whxcn a response xs desired 
should be entered xn these colunns xn the same x?ay as the 
frequency xncrement between prxntouts was entered., 

Txme response 

An xntegratxon step sxse may be entered xn these columns, wnose 
value governs the numerical integration producing the txme 
response. If omitted, a step size guaranteeing roughly six- 
figure accuracy of the resulting response will be used. If 
present, either this or the next smallest step size which 
evenly divides the print interval will be used. If the response 
to external sampled input is requested, the sampling interval 
of the external signal must appear m these columns, and the 
print interval must be a multaple (>2) of it. 

5.2 END CARD 

The last card in the input contains an asrensk in column one 

and has the free form format otherwise as 'indicated, 

* END. 
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CHAPTER VI 


EXAMPLES OP USE 

6.0 A -PULSE FORMING NETWORK 

A circuit consisting of two twin-tee networks connected by a 
negative impedance converter is mdacated by Figure 4. STICAP input 
cards which might be used to analyse this circuit are indicated 
below. The node numbering m the elements cards group agrees with 
that given m the figure. A zero- valued capacitor CZ is used to 
establish a branch at the output across which the output voltage is 
available. First employed by Pottle in the original COMAP program, 
this circuit is a realization by Antreich and Gleissner of an 
optimum pulse forming filter proposed by Jess and Schussler.^ The 
filter should have at most a one percent step response overshoot 
and a one percent stop band frequency response. 

6.1 GEAR MODE ANALYSIS? USEFCN OPTION 

\ 

The user's control cards needed to achieve an analysis of the 
previous circuit using the Gear integration routines are listed 
below m the order they would appear m the user's deck. In order 
to drive the network with the voltage source Vj_ as a sinusoidal 
input the following CUTOUT cards were used to alter USEFCN: 

CUTOUT 7700000 9700000 

U(l) *» SIN{5.0*T) 


3 K. Antreich and E. Gleissner, "Uber die Realisierung von Impulsfil- 
tem durch aktive RC-Uetzwerke , * Arch. Bisk. Ubertr ., Vol. 19, 1965 
pp. 309-316. 

4j. Jess and H„ W. Schussler, "On the Design of Pulse-Forming Net- 
works," IEEE Trans . , Vol. CT-12 (September 1965), pp. 393-399. 
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Before exercising this option the above sequence numbers should be 
verified by obtaining a listing of the program. (The CUTOUT option 
may be a feature peculiar only to the Langley computer? if so, other 
methods of altering USEFCN needs be employed.) 

The first card appearing below is the user's title card: 
ANTJREICH AND GLEISSNER - SCHUESSLER PULSE FORMING NETWORK 


^ELEMENTS 

R1 5 1 

R2 12 

R3 4 0 

R4 2 3 

R5 3 6 

RS 6 7 

R7 7 8 

R8 10 0 

R9 8 9 

Cl 10 

C2 2 0 

C3 3 0 

04 4 1 

C5 4 3 

C6 6 0 

C7 8 0 

C8 10 7 

C9 10 S 

CZ 9 0 

VI 5 0 

£D1 3 0 

*OUTPUTS 
VCZ 
*GEAR 

*SOURCE ORDERING 
VI 

*RUN CONTROLS 
STOP TIME =20.0 
OUTPUT DENSITY = 10 
*END 


1.0 

3.118 

3.121 

2.86 

0.9005 

3.127 

3.257 

0.6412 

11.21 

2.211 

0.749 

9.809E-2 

0.1713 

0.1867 

0.459 

0.196 

0.5972 

0.1735 

0.0 

-2.0 IR5 


Since the circuit is not stiff, this analysis could also be more 


efficiently achieved in the Gear mode by use of Adam's integration 
techniques. In this case the following card must be included in the 
RUN CONTROLS card group: 
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ADAMS INTEGRATION 

Other cards controlling the integration step size and output data 
might also be prepared. 

6.2 MATRIX MODE ANALYSIS; SOLUTION EQUATIONS PRINTED 

The analysis of the sinusoidal driven pulse forming circuit by 
matrix mode routines may be achieved by using the same title card 
and elements cards as m section 4.1, with the remainder of the 
user's input control cards replaced by the following: 

^OUTPUTS 

VCZ 

*MATRIX, SO 
*RUH CONTROLS 
STOP TIIIE = 20 
PLOT INCREMENT =0.5 
^SOURCE DEFINITIONS 
VI = SIIT (5 . 0*T) 

*END 

In this instance the closed form solution for VCZ as a function of 
time is also obtained. 

6.3 CORNAP MODE ANALYSIS; SAMPLED DATA INPUT 

A CORNAP mode analysis of the pulse forming network with a 
sinusoidal driver requires the same elements cards as in section 4.1; 
the remainder of these cards would be replaced by the following : 

*OUTPUTS 


VCZ 



40 


*CORHAP , STATE EQNS , TRANSFER FCIIS 
Control Card (See below for contents) 

DATA CARDS 
*END 

The contents of the nonblank columns of the control card mxght bes 

Columns Contents 


1 

T 

3-4 

VI 

10-12 

100 

28-31 

0.10 

33-35 

VCZ 

41 

X 

58-61 

0.01 


The data cards contain, 6 samples per card in (perhaps) F 12.4 
format, values of sm 5 t on the interval zero to ten seconds, at 
increments of 0.01 seconds. A smaller increment could be employed, 
in' the interests of accuracy of the numerical integration. A chief 
encumbrance of this mode of analysis is the preparation of numerous 
sampled data cards necessary for an accurate numerical integration. 
The options exercised on the mode card effect printing of the state 
equations, transfer functions, poles and zeroes of the network. 

6.4 S0I1E PROGRAM RESULTS 

In Figure 5 we exhibit for the pulse forming circuit with sinu- 
soidal driving function the voltage VCZ as a function of time, as 
obtained by execution of the STICAP program m the various modes, 
using the control cards of Sections 4.1, 4.2, 4.3. Figures 6 and 7 
contain the corresponding results for the step and impulse responses 
of the pulse forming network. Output from the three modes appears 
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to be mutually consistent* However, at values of time for which the 
output "voltage is nearly zero# the matrix mode printed results do 
not appear to have as many significant figures of accuracy as do the 
Gear mode results, suggesting loss of significant digits in this 
computational mode. The problem has not been thoroughly investigated 
to determine its significance. 


6.5 A STIFF CIRCUIT 

A stiff circuit which has caused problems m the first generation 
version of SCEPTRE^ is exhibited m Figure 8. The node numbering 
indicated agrees with that used xn preparing the control cards below. 
The system matrix for this network has the eigenvalues (poles of the 
system) -.5 + lOOi, -_Jl, -10. 

The elements card group and some possible outputs for the net- 

l 

work are: 

* ELEItENTS 

RS 1 2 1.0 

R1 2 3 1.0 

R2 2 5 1.0E + 4 

Cl 3 4 1.0E - 4 

C2 5 6 1.0E - 2 

LI 4 0 1.0 

L2 6 0 10.0 

VI 1 0 

* OUTPUTS 

VC1 

VC2 

IL1 

IL2 


6 C. H. Cooke and E. Young, "Numerical Integration and Other Techniques 
for Computer Aided Network Design Programming," NASA CR-111837 „ 
January, 1971. 
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The Gear and Matrix modes were used to obtain time response for 
this stiff circuit/ with the sinusoidal independent source 

VI = sm 5 t. 

Plots of the output voltage across capacitor Cl, VC1 versus time, 
appear in Figure 9. 

In Figure 9(b) the data is plotted on a more microscopic scale, 
utilising more data points, to illustrate the parasitic high fre- 
quency effects. These effects are most readily discemable m the 
first 0.5 seconds of time response. Tjhen such high frequency effects 
are of interest, many data points need be plotted, and the Gears in- 
tegration routine should be restricted with a stepsize UMAX small 
enough to give an output data density sufficient for discernment of 

i 

such effects. Whether or not a circuit is stiff may be detected by 
examining the poles of the network (eigenvalues of the system matrix) , 
obtainable by use of the transfer functions option, available in all 
modes. Large magnitude left halfplane poles indicate a stiff circuit. 

6.6 A CIRCUIT WITH SOURCE DERIVATIVES 

A circuit characterized by source derivatives in the state 
equations (due to a loop containing only inductors and capacitors) 
appears in Figure 10. The elements cards and possible output cards 
are? 

* ELEMENTS 

LI 1 2 

Cl 1 2 

C2 2 0 

VI 1 0 


1.0 

1.0 

3.0 



- 47 - 



- 0.1 
- 0.2 
- 0.3 
- 0 . 4 
- 0.5 
- 0.6 
-0 7 
- 0.8 
- 0.9 


- 1.0 









- 50 


* OUTPUTS 


VC1 

VC2 

ILl 


Figure IX depicts the output voltage VC1 resulting from the source 



Figure 10. A Circuit Characterised by 
Source Derivatives in the 
Circuit Equations 


6.7 OUTPUT LISTINGS 

This section contains a program output listing for the pulse 
forming circuit previously discussed. The listing illustrates the 
manner in which the matrices of the state equations and state output 
equations , transfer functions, gain constants, zeroes and poles of 
the network are printed. This listing was obtained using the CORNAP 
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mode vi th the sinusoidal driving function sm 5 t. The page contain 
mg the comment 

** The name VCZ does not appear on the output list ** 
is the result of misusing the CORNAP program. In specifying the 
program output , the voltage through CZ , the name VCZ was placed in 
columns 33-35 of the COKWAP control card (see page 40) . However, 
these card columns should have contained the name CZ, the same name 
as was used m specifying this capacitor m the elements description 
without the V appended when specifying the outputs by the outputs 
description card group. 
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ANTREICF AND GLEISSN6R - SCHOESSLER PuLSE FCRrtl Kb NETWORK -64. 10-10. S 
ELEMENTS 


R 1 

5 i 

1.0 

R2 

1 2 

3.118 

R3 

4 0 

3.121 

R4 

2 3 

2.86 

R5 

3 6 

C.9C05 

R6 

6 1 

3.127 

R7 

7 8 

3.257 

R 8 

10 C 

0.6412 

R9 

8 5 

11.21 

Cl 

1 C 

2.211 

C2 

2 0 

C.749 

C3 

3 0 

5.8Q9E-2 

C4 

4 1 

Ca 1713 

C5 

4 3 

0.1867 

C6 

6 C 

0.459 

C7 

6 Q 

0.196 

C8 

1C 7 

C. 59 72 

C9 

10 9 

C. 1735 

cz 

9 0 

0.0 

VI 

5 C 


IDX 

3 0 

-2-0 I 


QLTPLTS 

vcz 

THIS NETNCPK HAS 8EEN SCALED FCR CC MPliT AT ION 8Y TFE FOLLOWING FACTORS 
FRECDENCY 1-OOOCOUOE+OO RACIANS/SEC. IMPEDANCE I.CCCOCCOE+OO OHMS 



ANTREICH ANC GLEI5SNER 


5CFU £ S S L ER PLL5E FORMING NETWORK -6 4„ 10 » 10. S 


TF£ (SCALED) ENTRIES CF THE: A MATRIX ARE 


STATE 

VARIABLES 


STATE 

VARIABLES 



Cl 

C 2 

C3 

C4 

Cl 

-6. 8815 173E-01 

2. 1 5667 E9 fc-01 

1-6062439E-01 

- 1. C3 18 169E-01 

C 2 

4.28L9547E— 01 

-8. 5501637 £-01 

4.6682290E— 01 

0. 

C3 

-1.2195434E+00 

l.Sfc£446CE+GG 

4. 1 3603 12E +00 

-9.40721S7E-01 

C4 

- io 172 1256E+0G 

9. 13S 8 >4 15E-01 

2.0732080E+00 

-1,-33178466+00 

C6 

Oo 

0. 

2-4 193773E + 0O' 

-Oi *' 

C7 

0. 

0 • 

0. 

0. 

ca 

0. 

0. 

0. 

0. 

C9 

G. 

0. 

C. 

0. 


C 6 

Cl 

C8 

CS 

Cl 

-2 -344402 SE— 01 

c. 

- 0. 

C- 

C2 

0. 

0- 

0. 

0. 

C3 

-6-03676 92E+ GO 

C. 

0. 

0. 

C4 

-3-0259631E+0Q 

0. 

0. 

0. 

C6 

-3-0161£7SE + CC 

1.2131571 £-01 

-5.0480674E-01 

2.7312178E-02 

C7 

2.841Q159E-G1 

-1.6696C64E+0C 

-1.CC96212E+0G 

-3. 758836GE-0 1 

ca 

-3-87587766-01 

-3.31355 52 E -01 

-7.604893 7E-01 

4.1 145689E-02 

C9 

7-225527 1E-C2 

-4 . 2462 52 1E-C1 

1.4 162654E-01 

-4.S40C058E-01 


THE (SCALED) ENTRIES CF THE B MATRIX ARE 


STATE 

VARIABLES 


SCURCE 

VARIABLES 


Cl 

Vi 

4-429 18C 8E— 01 

C2 

C- 

C3 

2.111 1348E-0I 

C4 

- I.2C888C4E— Cl 

C6 

0 . 

C7 

0 . 

C8 

C. 

C9 

0. 


TFE (SCALED) ENTRIES OF THE C MATRIX ARE 



CUTPLT 

variables 



STATE 

VARIABLES 


V Cl 

Cl 

c. 

C 2 

C. 

C 3 

0. 

C4 

0 . 1 

V Cl 

C6 

la 40 S3 18 IE- 01 

C 7 

1.74123 

ca 

5SE-C1 2. 7M5442E-01 

C 9 

-9.6C79902E-01 


TEE ( SCALE C ) EMP1ES CF THE D MATRIX ARE 


CUT P LT 
VARIABLES 


SOURCE 

VARIABLES 
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ANTREICP ANC GL£ ISSNER - SCHOESSLER PULSE FCRMING KETHCRK -64. 10. LO. S 
TRANSFER FUNCTION CRITICAL FREQUENCIES ( SCALED! 


CUTPUT VARIABLE - V CZ 
SOLACE VARIABLE - VI 

GAIN CONSTANT IS 7. 1 7 7847CE-G2 

j 


PCLE 1 

POSIT IUNS 


ZERO 

POSITIONS 


REAL PART 

IMAGINARY PART 

CRCE R 

REAL PART 

IMAGINARY PART 

CRDE 

-4. 48 59Q94E-02 

9.1 549708 E-0 L 

1 

-6.0964382E-06 

1.29 45752E + Q0 

1 

—3*666881 4E-01 

6.2457201E-C1 

1 

- 8. 9 b 006 1 8 E-0 l 

0. 

1 

-8.95C0618E-G1 

0. 

1 

-2.0220162E + 00 

0 . 

1 

— 2. 02 20 16 IE -e-CO 

C. 

1 

l.y922384E-04 

1.0 19909 3 E+OQ 

1 

-4.9CE4538E-01 

2 *0 6103 36E—C I 

1 
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ANTREICH AND GLEISSNER - SCHUESSLER PULSE FORMING NETWORK -64. 10,10. S 


TRANSFER FUNCTION CRITICAL FREQUENCIES (SCALEO) 


CUTPUT VARIABLE - V CZ 
SOURCE VARIABLE - VI 


GAIN CONSTANT IS 


7.173847CE-C2 


POLE PCS IT IONS 
REAL PART IMAGINARY PART 


CRDER 


REAL 


ZERO POSITIONS 
PART IMAGINARY 


PART CRI 


-4.4E59C94E-02 

9.1 5497G8E-Q1 1 

-6.09643 82 E— 06 l . 29 9 5 752E+00 

-3.6668814E-C1 

6.2457201E-C1 1 

-8.95C0618E-01 0. 

-8.95C0618E-UI 

0. 1 

— 2.02 20 L 626+00 0. 

-2.0220 16IE+CC 
— 4.9Q54538E-Q1 

0. 1 

2.081G336E-C1 1 

1. 99223 84£ -04 1.0 199093E+00 

THIS NETWORK HAS BEEN SCALED FOR COMPUTATION BY THE FCLLCWING FACTORS 

FREQUENCY l . OOOGOOOE+OO RADIANS/SEC 

UNSCALEC TIRE 
INTEGRATION STEP 

. IMPEDANCE 1. QCCOOOCE+QQ OHMS 

RESPONSE 

SIZE 2.0000E-02 SEC 

external input 1 

SIGNAL SAMPLING INTERVAL I.G0G0E-02 SEC 

TIME 

INPUT 

CUTPUT (\jLZ) 

0 . 

0. 

0. 

i.OGCCE-Ql 

4.7943E-01 

5.6497E-G5 

2.0QCCE-01 

8.4147E-01 

A. 165QE-04 

3.00CCE-01 

9.9749E-C1 

1.2624E-03 

4.<!)000E-Ql 

9.093CE— 0 1 

2.6154E-03 

5.CCCCE-01 

5.9847 E-0 1 

4. 3398 E- 03 

6.GOOOE-01 

1.4112E-01 

6. 1835E-03 

7.GQCQE-01 

-3.5C7HE-C1 

7.8467E-03 

8.00CQE-01 

-7.568CE-01 

9.0615E-03 

9.0UCOE-01 

-9.77E2E-C1 

9.6622E-03 

I.OOCCE+CG 

-9.5892E-01 

9 .63006-03 

1.10CCE+00 

-7.055AE-CI 

9. 1015E-03 

1.2CCCE+00 

-2.7942E-C1 

8.3379E— 03 

1 .3QC C£ * C 0 

2.1512E-01 

7.6633E— 03 

1.40CCE+Q0 

6.56S9E-01 

i. 38 73E-C3 

I • 5QC CE+OQ 

9.38CCE— 01 

7. 7301E— 03 

1.60CCE+00 

9.893EE-0 1 

8.77C3E-03 

1.7QCCE+0Q 

7.9849 E-0 1 

1. 04256-02 

1. 30 C CE+OQ 

4.1212E-01 

1. 2472E-02 

I.9CCCE+00 

-7.5151 E-02 

1.46C1E-C2 
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~ 2* CGC CE + OC 
2 .lOCCE+uO 
2.2GCCE+QQ 
2 .30GCE+00 
2.4QCCE+CC 
2.5QC CE+OO 
2o6CCCE+0U 
2.7GCCE+QQ 
2.SOOQE+OU 
2.9QCCE+00 
3.GGC0b+GG 
3. IQCtE+OG 
3 *2GC CE + UO 
3.30CCE+G0 
3.40CCE+0G 
3 .50 C CE+UQ 
3.60CCE+D0 
3.7GCCE+00 
3.8GCCE+G0 
3.90QCE+00 
4.00QCfc+0Q 
4.1GCCE+U0 
4.200QE+00 
4.3CCCE+Q0 
4. 40 0 CE +QO 
4*50 C CE+OO 
4.6QCCE+0G 

4- 7000E+00 
4* 80 C CE +00 
4.90CCE+00 
5.CCC CE+UO 
5.10CCE+00 

5- 200CE+UQ 
5.30CCE+00 
5.4u C CE +0 0 
5 . 5CC CE + OO 

5 *60 C CE +00 
5.70CC6+00 
5.80CCE+0G 
5*90CCE+00 
6.00C CE+OO 
6.10C0E+00 
6.20CCE+00 
6.30CCE+00 

6 *40C CE + OO 
6.5JCCE+00 
6 *600 CE +00 
6.70CCE+00 
6* 80C CE +00 
6*90CCE+00 
7.00CCE+00 


-5.44U2E-01 
-8.797CE-01 
-9.5595E-C1 
-8*7545 E— 01 
-5,36575-01 
-6.6322E-C2 
4.2G17E-01 
8.0378E-01 
9.906 IE— 0 1 
9.349CE-01 
6.5C2SE-01 
2.064 7E-0 1 
-2 . 879CE-0 1 
-7*1179 £-0 1 
-9.614CE-01 
-9.7563E-0 l 
-7.5C55E-01 

— 3 .424 6E-0 l 
1.4988E-01 
6. C554E-0 1 
9.12S5E-01 
9. 9683E-C 1 
8. 3666E-0 1 
**. 7 164E-0 1 

-8.85 13E-C3 

— 4 .8717 E— 0 1 
-8. 4622E-C 1 
-9.58C6E-01 
-9. 055 E E-0 1 
—5.91 36E-0 1 

— 1*32 3 5 E-0 1 
3.59U6E-C1 
7. 6256 E- Cl 
9 • 7S36E-0 1 
9. 5638 E- 01 
6.9524E-01 
2. 709 IE— 0 1 

-2.2376E-01 
-6.6362E-0 1 
-9.4103E-G1 
-9*bbC3E-0 1 
-7.9313E-01 
-4 * 04U4 E-0 1 
8.3974E-G2 
5.5143 E-0 1 
8.8387E-01 
5.499 1E-C1 
ti.7114E-01 
5.29C8E-01 
5.7467E-02 
-4.2818E-C1 


1.6494E-02 
1 * 7898E-02 
1 . 8690E-02 
1.8902E-02 
1 <>87 16E-02 
1* 84 19E-02 
1. 83 26E-Q2 
1*87 C9E-02 
1.9 723E-02 
2.1371E-02 
2. 3500E-C2 
2.5837E-02 
2.8057E-C2 
2.9859E-C2 
3* i C39E-02 
3* 1 539E-02 
3.I461E-C2 
3. 1040E-C2 
3*0 586E- 02 
3* 04C7E-C2 
3* 073 2E-02 
3.1656E-C2 
3* 31 1 56-C2 
3*48996-02 
3.67C8E-02 
3. 8219E-02 
3.9168E-02 
3.9414E-C2 
3.8973E-C2 
3. 8013E-02 
3 . 68 12E-G2 
3* 5694E-02 
3 . 4944E-02 
3 .47 4 IE-02 
3*511 5E-C2 
3.59386-02 
3. 6957E-02 
3.7854E-02 
3* 8326E-C2 
3*8159 E- 02 
3. 2282E-C2 
3.5/8 1E-C2 
3.3682E-C2 
3. 1896E-02 
3.014 1E-U2 
2. 686 7E-02 
2.81 95E-C2 
2.8C86E-G2 
2 * 8356E-C2 
2.8716E-C2 
2 . 8 846E-02 
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7.1CCCE+QU 

7.20CGE+00 

7.30CCE+00 

7.40CCE+00 

7.50UCE+00 

7.60CC6+00 

7.7GCCE+OQ 

7.80CCE+00 

7*90CCE+00 

8.GGCCE+00 

8* ICC CE + OG 

8 .200GE+GC 

8.3CCCE+00 

8 .4QCCE+0Q 

8.500 GE+OU 

8.60C Cfc + OQ 

8.7QLCE+QG 

8.8GCCE+00 

8.90CCE+0G 

9-OOCCE+OO 

9.1QCCE+0Q 

9.20CCE +00 

9.30CCE + GC 

9.40CGE+00 

9.5GC CE + OO 

S.6GCCE+0G 

9.70GCE+0Q 

S.8OCCE+G0 

9.9ULGE+0J 

l.CGCCE + O l 


-8.0902 E--0 1 
- 9. 91 76E-Q 1 
-9. 3172 E- 01 
-t>.4354E-Cl 
- 1 . 9 78CE-Q 1 
2.9637E-01 
7.17S7E-01 
9.638CE-G1 
9.736 5 E-0 1 
7.451 IE-01 
3.3415E— Cl 
-1.5862 E-0 1 
-6. 1256E-C 1 
-9. 1652 E-G 1 
-9.96C9E-Q1 
-8.3177E-01 
-4.6382E-Q1 
1 .7702E-Q2 
4.5489E-G1 
8.509CE-0 l 
9.585SE-C1 
9. 0179 E- Cl 
5. 842 C E-0 1 
1.2357E-C1 
-3.673 IE-01 
-7.6825E-G1 
-9.8111E-01 
-9.5375E-C 1 
-0.92t€E-Cl 
-2.6237E-CI 


2.8478E-C2 
2.7455E-C2 
2.5777E-C2 
2.3597E-02 
2.116 7E-C2 
1.88 74E-02 
I. 69 55E-C2 
1.5632E-Q2 
I. 4958E-C2 
1. 4829E-02 
1.5GG7E-02 
1. 5182E-C2 
1. 5047E-02 
1*437 5£— 02 
1. 30 73E-02 
1.1209E-02 
£. 596QE-03 
6. 736QE-Q3 
4* 7517E-03 
3.3060E-G3 
2.5385E-G3 
2.43 18E-C3 
2.8164E-G3 
3. 4 125E-Q3 
3. 899 IE-03 
3. 9925E- 03 
3.5 1 67E-03 
2.4458E-G3 
9.1146E-C4 
-8.3001E-04 
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The next listing depicts Gear mode output for the pulse forming 
circuit when step response is requested. Observe that the input run 
controls cards used are listed as part of the output; as well as the 
elements and outputs card groups : 
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ANTREICH AND GLEISSNER - SCHUESSLEP PULSE FORMING NETWORK - 64,10*10.$ 
ELEMENTS 


Rl 

5 

l 

1-0 

R2 

1 

2 

3. 1 18 

R3 

4 

C 

3.121 

Rf 

2 

3 

2.E6 

R5 

3 

6 

C. 9005 

R6 

6 

7 

3.127 

R7 

7 

8 

3.257 

R 8 

10 

0 

0.6412 

R9 

0 

9 

11.21 

Ci 

1 

0 

2.21! 

C2 

2 

0 

0.749 

C3 

3 

0 

9.F09E-? 

C4 

4 

1 

C »' 71 3 

C5 

4 

3 

G. 1 867 

C6 

6 

0 

0.459 

C7 

8 

0 

0.* 96 

cs 

10 

7 

0.5972 

C9 

10 

9 

0." 735 

cz 

9 

0 

O.C 

VI 

101 

5 

3 

0 

0 

-2.C IR5 


OUTPUTS 
VC 7 

THIS NETWORK HAS BEEN SCALEO FOR COMPUTATION BY THE FOLLOWING FACTORS 


FREQUENCY 1 .OOCOOOOE+OO RADIANS/SEC 


IMPECANCE l .OOCOOCOE+OO OHMS 
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ANTREICH AND GLEISSNER - SCHUESSLER PULSE FORMING NETWORK - 64. 10. 10. S 


THE (SCALED) ENTRIES QF THE A MATRIX ARE 


STATE 

VARIABLES 


STATE 

VARIABLES 



Cl 

C2 

C2 

C4 

Cl 

-6.381 51 75E-C1 

2.1586789E-01 

1.6062439E-01 

-1 .03181696-0 1 

c^ 

4.28I954-'E-C1 

-8.3501837E-01 

4.668223CE-C1 

C. 

C3 

-1.21354 34 E+ 00 

1.9684460E+00 

4.136C312E+C0 

-9.4072197E-01 

C4 

-1 .1721256E+00 

9.1338415E-01 

2.C73298CE+00 

-1'. 3317 846E + 00 

C6 

0 . 

0. 

2.4193773E+00 

C. 

C T 

0. 

0. 

0. 

0. 

C9 

0. 

0 . 

0. 

c. 

C9 

0 . 

0. 

0. 

0 . 


C6 

C7 

C8 

C9 

Cl 

-2.3444C23E-C1 

0. 

0. 

0. i 

C2 

0 . 

0. 

0. 

0. 

C3 

-6<. 0367632 E+CC 

0. 

0. 

0. 

C4 

-3 .0253631 E+CO 

0. 

0. 

0 . 

C6 

-3 .01 61875F+00 

1 . 21 31 571E-01 

-5.048C674E-01 

2.7312178E-02 

C 7 

2.84’ 0153P-01 

— 1 . 66960 6 4E +00 

-1.C096212E+00 

-3 .7588360 E— 01 

C3 

-3.87 C 8776E-C1 

-3. 31355926-01 

-7.6048937E-01 

4 • 1145689E-02 

C9 

7.2255271 E-C2 

-4. 246293 IE-01 

1 .41 62654E-C1 

-4.9400G58E-01 


THE (SCALED) ENTRIES OF THE B MATRIX APE 

STATE SOURCE 

VARIABLES VARIABLES 


Cl 

VI 

4.42318C8E-01 

C 2 

0 . 

C3 

2 .11 1 1348E-01 

C4 

-1 « 2C 688C4E— Cl 

C6 

0 . 

C7 

0 . 

C8 

0 . 

C9 

0 . 


THE (SCALED) ENTRIES OF THE C MATRIX ARE 



62 


OUTPUT STATE 

VARIABLES VARIABLES 




Cl 

C2 

C3 

C4 

l v 

CZ j 

0 . 

0 . 

0 . 

0 . 


1 

C 6 

C7 

C8 

C9 

V 

CZ 

1.4C531 81E-C1 

1.7412359E-01 

2.7545442E-C1 

-9.6079902E 


THE (SCALED! ENTRIES OF THE D MATRIX ARE 

OUTPUT SOURCE 

VARIABLES VARIABLES 


V CZ 


0 


VI 



63 


ANTREICH A NO GLEISSNER - SCHUESSLER PULSE FORMING NEThCPK - 64.10. 10. S 


TRANSFER FUNCTION CRITICAL FREQUENCIES (SCALED) 


OUTPUT VARIABLE -VC l 
SOURCE VARIABLE - VI 


GAIN CONSTANT 

IS 7.177 8470E 

-02 




PCLE 

PCSITICNS 


ZERO 

POSITIONS 


REAL PART 

IMAGINARY PART 

ORDER 

REAL PART 

IMAGINARY PART 

CRDE 

-4.4859094E-02 

9 . 15497C 8E-01 

1 

-6.0564382E-06 

1.2945752E+00 

1 

-3. 66 6 881 4 £-01 

6.2457201E-01 

1 

-8.9500618E-C1 

0 . 

1 

-8.0,50061 8E- 01 

0 . 

1 

-2.0220162E+OC 

0 . 

1 

-2 . 0220161 E+00 

0 . 

1 

1. 99223846-04 

1.C1990S3E+00 

1 

-4.QQ54538E-01 

2 . C 8 10 336E-01 

1 
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RUN CCNTPCLS 

EPS = 1 . 30CCCE-C4 

PRINT DENSITY = 3 

INITIAL TIME = 0. 

PRINT STARTING TIME = 0. 

HINIT = l.CCOOOE-04 
HMIN * 0. 

STIFF INTEGRATION 
HMAX = 2. CCCCOE+OO 

STOP TIME = 2.0CC00E+01 

STEP RESPONSE 
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ANTRgrCH AND 
TIME 

0. 

1.161C4E-02 
5. 29586E-C2 
1 .08541E-01 
1 .63856E-01 
2,01?3TE“Ql 
2.32724E-01 
2»°4326E“Q1 
3.71274E-01 

H. 39755E-01 
5.21930E-01 
6.D7320E-01 
e> . 97CC8E-01 
7.90698E-QI 
8.91607E-01 
1 » 00403 E-*-0D 

I. 125766+CO 
U25278E + 00 
1 .382C7E+00 
1.51220E+00 
1 . 64242E+00 
1.77235E+0 0 
I.90I75E+00 
2 . 03054E + 00 
2. 15877E+00 
2.28659E+00 
2 . 41 874E *00 
2 , 674S7E+00 
2.9895CE+C0 
3.38514E+00 
* o S8472E*C0 
4.492588*00 
5.15459E+00 
5.90588E+00 
6.6571 8E + 03 
7.40847E+00 
8.1 5977E*-00 

8.91 l 068*00 
9.66236E*00 
1.041 37E*0l 


GL8ISSNER - 5CHUESSLER PULSE FORMING NETWORK - 64.10. 

VCZ 

0 . 

4.77371E-C6 
9.73736E-05 
3.96535E-04 
8.75283E-04 
1.3C866E-03 
1 »7C8£ 4 £-03 
2.63862E-03 
4.025C3E-03 
S.45897E-03 
7.39103E-03 
9.62655E-03 
1.22051E-02 
1.51354E-02 
1 .6552CE-02 
2.26723E-02 
2.75121E-02 
2.3CC28E-02 
3.908C3E-02 
4.57345E-02 
5.29751E-02 
6.08256E-C2 
6.93056E'02 
7.9455CE-C2 
3.829398-02 
9.88619E-02 
1.1C616E-01 
1 . 358 86 E— 01 
3 . 7 J482E-C2 
2.23322E-01 
2.9CC58E-01 
4 .032978-01 
5.25C62E-01 
6 .622 36E-01 
7.9567CE-C1 
E.8492CE-C1 
9.542 62 E-Ol 
8. 93620 E- 07 
1 .008 16F*00 
1 .0C656E+-C0 
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ANT RE ICH AND GLFISSNER - SCHUESSLER PULSE FORMING NETWORK - 64. 10. 10. S 

TIME VCZ 

1.T1650E+C1 9. <383946-01 

l .19850E+01 9 . 91 1C0E-01 

1 . 28393E+01 9.906C9E-01 

1.38611E+01 9.9 8443E-01 

1.49667E+01 1.0C782E+0C 

1 .607226*01 1 • 0C871E+OO 

1.71 777E+01 1.0CC88E+00 

1 .82832E+01 9 .925 63E-01 

• 93887E+01 9. 91795 E- 01 

• OOOOOE+Ol 9.95055E-01 


So I y ioh Co m pie fc 
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08/17/71 LRC f N SCOPE 3.0 64C0Z-131K 06/02/71C 

15.08 . 16 .TER 1 31 7 . 

15.08.17. JCB,i»2CC, 700CC. A3031 R 

15.08. 17. HB105 BLOG 12C2 CENT 

15.08 .17. USER. BAVUSO, SALVATORE J 0000 

15.08.17 .44430 15010 

15.08.17. NCM AP. 

15. 08. 17. FETCH (A 3^ , SPR ZOl , BI NARY) 

15. 08. 23. TIME BG ATTACH 

15. 09. 09. TIME ED ATTACH 

15. 09. 28. END FETCH 

15.09.28. REWINCICUTPLT ) 

15.09.30.BNF ILE . 

15. 10. 34. STOP ' 

15. 10. 35. CPU 7.134427 SEC. 

15.10.35. PPU 1 16.859840 SEC. 

15. 10. 35. TL. *77 SEC. 

15. 10. 35. DATE C8/17/71 

16.07.15. TER13 1 7. 215 LINES PRINTED. LP26 
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The next listing exhibits typical matrix mode output for the 
pulse forming circuit with the sm 5 t driver. Observe the form of 
the solution equation giving VCZ as a function of tine. The notation 
1. 31040E-09EXP (-2 . 02202E + 00T) + ... etc. 
is to be interpreted as 

1.3104 X 10" 9 e” 2 * 022027 

where T is the time variable. The T appearing in sm and cosine 
arguments is also to be interpreted as the time variable, and not a 
portion of any floating point exponent which may precede it. 




ANTREICH AMI) GLEISSNER 


ELEMENTS 


1.0 


3,118 


R3 4 U 3.121 


2 . So 


R5 6 o.so 05 


3.1?7 


NQ.......13JL9. 

SCHUESSLER PULSE 


FORM I MG NFTUORK -64.10.10 



8 

10 7 


lu 9 


Cl 9 0 
“ 7 I~ 5“3 


ini iu 


OUTPUTS 


VC L 

TH I o HE n<Ur\K 


0.3072 


0.1735 


0. 0 


- 2.0 


HAS 3CCM SCALED FOR COMPUTATION 


T i ! r FOLLOl INS FACTOR 


FREQUEN CY 1 . UUU UU OUF + 0 U RAD I AM 5/SCC . I <PE OAK f 1 . QOu JO QOF+QO O’M 
AN f REICH AID *Tle' 1 sS"E jT - 'SCMtJ ‘-‘SS LER _P J LSE F 5® 7|V NE TWORK -64.10.10. 

RUil CUiiTUoL ■> " “ _ 

IN IT! \L Tl.iF =«)."” 

PR I .IT STARTf'n Tl i! = ~ ' 

iHli 1 E < ) L ~ POINTS PRliT-0 ■_ 4 o 

PRINT LITE iVAl '= ' 5.i)u‘i)0"t)F-ul 
STOP T I It = 2.J^OuJr + U3' 
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AU fN lTlAL CONDITIONS ZERO 

SOU KG E DEFIN ITI ONS ' "Z I”_" . ", ‘ ~ ” ~ 

' S j > I EL Ful'ICTt ON HAGtrrfiJDF " I.ToIjWe+Uu AN 0“ AN OU LAL VELOCITY 5.U0"CP 
ANtKElCH AND GLEISqNER - SCH UESSI.FR JP ‘LSE__FOP H M Q HFTMORK -64.10.1U.S 
- THE FOLLOWING AkE THE CLOSED " FUR. j SOLUTIONS' 


VCX 

"2. o 6 19 IF -0 3 S 1 i K 5.G 0 ()u U_E+UU T + 3 3 2 85F- 3 U 

~~ ~+ II 3 1 Q 4 0 E - U J EX P C -T", 02 2T2~ E + Q u" T ) - 8 5-7)1?) E - 0 IT ) 


+ 4, 

37b5 4E-01EX P( -4 .UU5 45E-0 IT) *S 1 U ( 2 .{»« 103E-01T + 6. 14925F-02 ) 

-i. 

ibl2uE-01EXP(-3.oh6b8r-UlT)-*-Sri( 6.2457 7F-G1T + 2.44476r-01) 

+ 4. 

36887E-0 3EXP(“4.4J>byiF-U?T)*S! , ]( 9 . 15497E-31T + 1.5D343E+00) 

- ANT.tEICH and 

jLE 1 _S6NFi\ - SCHUISSLEk PULSE F.iR.UKO NFTWORK -64.10.10.5 

T ! i IE 

VC/. 

0. 

2 . &3b7SE-lb 

5. JOOQOE-Ol 

“7 . 31503E-04 

1. UOQOGF+UU 

2.11766E-03 

1.50OUUE+U0 

-2 . 6bl7 9E-03 

2.0u0U0E+00 

2. 1473UE-03 

- 2.5UO00F+UU 

-7.78797E-04 

5 . uOuuuE+UO 

-0. J9442F-J4 

3.50UUOE+UO 

2.2l9’JbE-J3 

4.uu0UUE+0u 

-2.65757E-05 

4.SUU00F+-0U 

2 .U3823E-U^ 

5 .Uu0UbF+ub 

-l. .U8264E-U4 

5.5(JOUUE + Ou 

. " 1 9 Go 3p 2E-03 APTr'TNT (ST 1 T 6 flTT* w 

b.oGOOOF+UU 

2 o 3124 jE-u 3 UixllallNiaULi if Avri!* jjj 

b . 5 UUUQF +Ul) 


7 . OUUUOE+ uU 

i. J2U20E-U3 

7.50uGbE+OU 

-4 .35U53E-04 

& . OOOUUE+Gi) 

-1. '2312E-05 

8 « 50 UGOE+OO 

2 „3d4o4E-0 5 

y.OooooE+ou 

-2 o bl410E-U3 

J.50UUUE-5-UD 

1.7U37UE-U3 


-2 . 5 j J2 bE -<J4 

WWIITI'IU'M 

-1 „ 37 723F -u3 

1. lUUUUE+Ul 

2 .4boo4F-J3 

1. 15UU0E+O 1 

-2 , 57504E-O3 

1.20UU0E+D1 

] . i>5 'J31E - J 3 

1. 250UUE+U T 

-8.3b54'*E-05 

1.5UUGUF+01 

-1.5 2b27 E-U 3 

1. 350UGF+01 

2 . 5 2 7 5 rf E - 0 3 

1.4G0UUE+')1 

-2 .524 j3E- J 3 

1.45UUUF+D1 

I. 517'>lF-o3 

1.50GOGE+01 

y . 29o4 JE-U5 

1.55UUOE+01 

-1 . LibobOf — U 3 

l.bUUUUE+01 

2.5 77 5UE-03 

1.G5U00E+01 

-? .4o5"ilE-U3 
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1.7U0 UUE+U1 1, 36 l i23E-l)3 __ 

T773qouf+ 0T "2Vb9T15E-J5' 
TTTiTDMF+tri "=TTS'UaFgr=D3 
1.8500UE + 01 2 * bl585E -0 3 

"1.9UU00ETJ1 — -2.39075E-03 

l.d5»0UF»Ul ~ "1.21432E- 0 3 


10/26/71 LUC On SCU RF 3. u 040 0 Z-151K 0 i/Ul/71 

15.au. 44. TFR1 319. 

15.10. 44. JUrf.l.jQO, 700U0. £3031 _R 

15.10, 44 .HR 10 5 BLO G 1 202” 

15 . 10.44. JSErt.OAVUSO, SALVATORE J ” 00 UP 

15.10.44.44 450 15010 

15 . lU.44.FETCTTA3031,bPRAU5, SOURCE ) 

15.1U.48.T1mE 0Q_ ATTAC H " 

15. 11. 29. Tl.iE _fj> ’ATT A CJL _ 

15. 15. u5. END FETCH 

15.15. Qo.riKSP,OUTPUT # 2. 

15.15. J7.dO.IAP. 

15.15,07. RUN ( S. , , SCF |_LE^ !) Jm) 

15 . 18 . 3b . SET in J h . _ _ 

15,18,57 . LGO . 

■ 15 . 19 . 14 .STv/P _ 

I 15.1-J.15. 0001005 0/S CALLb 

15. 10. 13. CPU 94 .0998 24 SEC , 

15, 19. 15 .PPU 142.5 4489b SE C. 

15. 19. 15. TL. 237 SEC. 


Ssji,,. 









APPENDIX I 


CIRCUIT SCALING 


Some basic concepts of the STICAP circuit scaling option tm.ll 
here be discussed. For further clarification one may consult Appendix 
C of the reference: 

William D. Stanley, Transform Circuit Analysis for Engineering 
and Technology , Prentice-Hall, 1968. 

The linearity of the circuits processable by STICAP allows the 
state-input and state-output equations to be algebraicly rearranged 
m the partitioned forms* 


a 

at 








V 

“V 


— 1 -1 

A 2 , B 2 

•m <■* <* 

+ 

f2 ! d 2. 

. H 


[V 


- i * 

]_ B 4 

1 ” . 

+ 

1 

' 

C 4 f 

U ' — 



(i) 


* 


( 2 ) 


Here the A, B, C, D subscripted quantities are constant matrices; and 
j is a time varying vector of state variables, 

Dj j is a time varying vector of the independent voltage sources 
or current sources , 

j is a time varying vector of user requested outputs. 
Moreover, the "v* and ,? i" subscripts indicate respectively voltages 
and currents of the network. 


i 

*For simplicity we consider only the case an which no source denva' 
tives are present. This renegade circumstance may similarly be 
treated . 
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Frequency Scaling 

User input of a frequency scaling factor 

FREQUENCY « w, oj f 1 

effects on the time scale the change of variable 

t — Xt * , X = o)~ x , 

where t* is the scaled time variable. 

The effect on the state-input equations is that state variables 
and independent sources are now measured in terms of the variable t* , 
and the A, B, C, D matrices m the scaled quantities are now the old 
A, B, C, D matrices multiplied by X. The form of the state-output 
equations does not change; however, the outputs are now expressed *.n 
terms of the scaled time variable t ! . 

As far as the user is concerned, the implications of using a 
frequency scaling factor are the followings 

(1) Sampled data inputs and Fortran defined or otherwise 

user supplied input data must be esqpressed m terms of 
the scaled time variable t'. 

(ii) Calculated time response outputs for which the unscale 
option is not (or cannot be) exercised are expressed m 
terms of the scaled time parameter t’ . 

(111) User controls on the numerical integration routines must 
be specified m terms of t’; i.e. stop time, print in- 
crement, start time, etc. 

(iv) Solution equations computed m the matrix mode x;ill be 
expressed m terms of the time variable t 1 . 
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Impedance Scaling 

Assume an input impedance scaling factor 

HJPEDANCE » kz, kz }* 1 

'is user supplied by means of the scaling option. In this instance 
the solution of equations (1) and (2) , but with the B and C subscript? 
matrices of the state-input ancl state-output equations replaced as 

i 

indicated below, is obtained. 

S( } _ -> jL B ( J 

(3) 

C { } } Jcz C ( J 

The implications of impedance scaling as seen from the new form of 
the state and output equations is that scaled outputs and unsealed 
outputs are not simply related (linearly) except m the following 
cases : 

(a) Only voltage sources 

(l) Scaled time domain voltage outputs from STICAP are 
true outputs. 

(n) Scaled current outputs must be divided by kz to 
obtain the true output. 

(in) Initial conditions input to the program must, of 

course, be scaled, inversely to that scaling of (i) , 

(n) . 

(b) Only current sources 

(i) Scaled time domain current outputs are true outputs, 

(n) Voltage outputs must be multiplied by kz to obtain 

true outputs. 



(ill) Initial conditions input to the program must be 
scaled, inversely to that of (i) and ( 11 ) . 

The case (a) corresponds mathematically to making on equations 
(1) and (2) the changes of variable 

(X 1 ) scaled = kz X x 

(Y . ) scaled = kz Y, 

-1* -i, 

with no changes of variable on Xy, Y v „ In case (b) the corresponding 
changes of variable are 

<X V > scaled - X v 

(Y v ) scaled - 5„. 

For the mixed case Stanley suggests the combining of voltage and 
current sources to obtain sources all of one kind. If this is not 
considered feasible the alternate below is advocated. 

For the mixed case the scaled equations produced by STICAP can 
alternately be obtained from the unsealed equations by either of the 
following changes of variable: 

I. (Xy) scaled = x v II. (X x ) scaled « kz X x 

(Y v ) scaled = Y v (Y x ) scaled = kzY x 

(U v ) scaled = kz - U v (U x ) scaled = kzU x 

Thus the user may use the scaling option for the mixed case, provided 

he scales {one but not both of) his current or voltage sources input 
to the program, as well as the corresponding initial conditions. The 
outputs are then to be interpreted by referral to (one but not both 
of) the proper output scaling equations {see case I and II) . 



